Strang splitting schemes for N-level Bloch models

In this paper, we define a splitting scheme for the [Formula: see text]-level Bloch model which makes use of exact numerical solutions sub-equations. These involve matrix exponentials want to avoid calculate at each time-step. The resulting is nonstandard and preserves qualitative properties equations. We explore compare numerically multiple ways implement it in particular take into account specific structure